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19 Abstract 


20 The objective of this study is to quantitatively evaluate the successive Tropical Rainfall 

21 Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) products and 


22 


23 


24 


further to explore the improvements and error propagation of the latest 3B42V7 
algorithm relative to its predecessor 3B42V6 using the Coupled Routing and Excess 
Storage (CREST) hydrologic model in the mountainous Wangchu Basin of Bhutan. First, 

25 the comparison to a decade-long (2001-2010) daily rain gauge dataset reveals that: 1) 

26 3B42V7 generally improves upon 3B42V6’s underestimation both for the whole basin 

27 (bias from -41.15% to -8.38%) and for a 0.25°x0.25° grid cell with high-density gauges 

. <*. ,0, % ,o 0.0 4%) , ^ pL t e— „ — 

29 coefficients (CC) (from 0.36 to 0.40 for basin-wide and from 0.37 to 0.41 for grid); and 2) 

30 3B42V7 also improves its occurrence frequency across the rain intensity spectrum. Using 

31 the CREST model that has been calibrated with rain gauge inputs, the 3B42V6-based 

32 simulation shows limited hydrologic prediction NSCE skill (0.23 in daily scale and 0.25 

33 in monthly scale) while 3B42V7 performs fairly well (0.66 in daily scale and 0.77 in 


34 monthly scale), a comparable skill score with the gauge rainfall simulations. After 

35 recalibrating the model with the respective TMPA data, significant improvements are 

r 

36 observed for 3B42V6 across all categories, but not as much enhancement for the already- 


37 well-performing 3B42V7 except for a reduction in bias (from -26.98% to -4.81%). hi 

38 summary, the latest 3B42V7 algorithm reveals a significant upgrade from 3B42V6 both 

39 in precipitation accuracy (i.e., correcting the underestimation) thus improving its 


2 




ACCEPTED MANUSCRIPT 


40 

41 

42 

43 

44 

45 

46 

47 


potential hydrological utility. Forcing the model with 3B42V7 rainfall yields comparable 
skill scores with in-situ gauges even without recalibration of the hydrological model by 
the satellite precipitation, a compensating approach often used but not favored by the 


hydrology community, particularly in ungauged basins. 


Keywords: CREST Model; A-Priori Parameter Estimation; H 
Evaluation; Precipitation Estimation 
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1 Introduction 


Precipitation is among the most important forcing data for hydrological models. It has 
been arguably nearly impossible for hydrologists to simulate the water cycles over 
regions with no or sparse precipitation gauge networks, especially over complex terrain 
or remote areas. Recently, the satellite precipitation products such as TMPA (Huffman et 
al„ 2007), CMORPH (Joyce et al„ 2004), PERSIANN (Sorooshian et ah, 2000) and 
PERSIANN-CCS (Hong et ah, 2004) are starting to provide alternatives for estimating 
rainfall data and also pose new challenges for hydrologists in understanding and applying 


the remotely-sensed information. 

The Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation 
Analysis (TMPA), developed by the National Aeronautics and Space Administration 
(NASA) Goddard Space Flight Center (GSFC), provides a calibration-based sequential 
scheme for combining precipitation estimates from multiple satellites, as well as monthly 


gauge analyses where feasible, at fine spatial and temporal scales (0.25°x0.25° and 3 
hourly) over 50°N-50°S (Huffman et al., 2007). TMPA is computed for two products: 
near-real-time version (TMPA 3B42RT, hereafter referred to as 3B42RT) and post-real- 
time research version (TMPA 3B42 V6, hereafter referred to as 3B42V6). 3B42V6 has 
been widely used in hydrological applications (Bitew and Gebremichael, 2011; Bitew et 
al., 2011; Khan et al., 2011a; Khan et al., 2011b; Li et al., 2012; Stisen and Sandholt, 
2010; Su et al., 2008), however, its computation ended June 30 th 2011 and 3B42V6 was 
replaced by the new version (TMPA 3B42 V7, hereafter referred to as 3B42V7), which 
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69 has been reprocessed and available from 1998 to present. Previously, 3B42V6 has been 

70 validated by several studies (Bitew and Gebremichael, 2011; Bitew et al., 2011; 

71 Chokngamwong and Chiu, 2008; Islam and Uyeda, 2007; Jamandre and Narisma; Jiang 

72 et ah, 2012; Li et ah, 2012; Mishra et ah, 2010; Stisen and Sandholt, 2010; Su et al., 2008; 

73 Yong et al., 2012; Yong et al, 2010), while the newly available 3B42V7 is evaluated in 

74 tropical cyclone systems (Chen et ah, 2013b) and the United States (Chen et al. 2013a 

75 and Kirstetter et al. 2013), it has not been extensively statistically and hydrologically 

76 validated in mountainous South Asian regions.. 

77 Therefore, the objectives of this study are designed (1) to evaluate the widely used 

78 globally-available, high-resolution TMPA satellite precipitation products over the 

79 mountainous medium-sized Wangchu basin (3550 km ) in Bhutan, and more importantly 

80 (2) to assess improve me ntssof the latest upgrade version (3B42V7) relative to its 

81 predecessor in terms of statistical perfonnance and hydrologic utility. Additionally, this 

82 study aims to shed light on the suitability of recalibrating a hydrological model with the 

83 remotely-sensed rainfall information. The remainder of this paper is organized as follows: 

84 Section 2 introduces the study area, the datasets used, and the methodology, including a 

85 brief description of the CREST distributed hydrological model and its upgrade to the new 

86 version (CREST Version 2.0). The results are discussed in Section 3, and then Section 4 

87 draws the conclusions of this study. 

88 
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89 2 Study Area, Data and Methodology 


90 2.1 Study Area 


91 The Wangchu Basin, with a total drainage area of approximately 3550 km is located 

92 within 89°6'- 89°46’E and 27°6'-27°5TN in the west of Bhutan (Figure 1). Wangchu 


93 Basin is the most populous part of the country with about 3/5 of the population living in 
1/5 of the basin area. The basin is equipped with one streamflow gauge at the outlet 
Chhukha Dam Hydrological station and five rain gauge stations. The soil types are 
dominated by Sandy Clay Loam (75.1%) and Loam (24.9%) based on the Harmoni z ed 
World Soil Database (HWSD vl.l) (FAO/IIASA/ISRIC/ISSCAS/JRC, 2009). The 

_ ...... ,« «. gJL .r — — r.,» 

99 (48.1%), woodland (17.8%), open shrubland (9.7%), wooded grassland (8.2%), grassland 

100 (7.6%) and other land-use types Hess than 10%) (Hansen et al., 2000). 


94 


95 


96 


97 


98 


101 


The northern periphery of the Wangchu Basin in the Himalayas has elevations over 


102 6000 m and maintains an annual snowpack. Lower portions of the basin are drastically 


. 


103 different and are subject to a summer monsoon from May to October (Bookhagen and 

104 Burbank, 2010). On average, the annual month with the greatest precipitation is July or 

105 August with 161 to 546 mm/month based on the five rain gauge station data shown in 


5 fi Aus 

5 (Tal 


3 1 

106 (Table 1), and the largest resulting streamflow occurs in June or August with 251m s' . It 

107 is possible that snowmelt contributes to a portion of this peak streamflow, but the 

108 majority is driven by the summer monsoon rains. In this study, neither the precipitation 

109 products nor the model explicitly deal with frozen precipitation. These are subjects 
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110 requiring additional investigation, especially in light of the forthcoming Global 


111 


112 


113 


114 


115 


116 


117 


118 


119 


120 


121 


122 


Precipitation Measurement Mission (GPM), which aims to quantitatively estimate frozen 
precipitation amounts. 

Insert Figure 1 about Here 
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2.2 In-situ and Satellite Precipitation Datasets 
2.2.1 Gauged Precipitation and Discharge Data 

Daily observed precipitation data are obtained from the Hydro-Met Services 
Department of Bhutan from 2001 to 2010 for the 5 rain gauge stations located within the 
Wangchu basin. In winter, frozen preetp.tation is reported in the form of water equivalent 
and computed by melting the ice/snow with hot water in the standard vessel and 
deducting the hot water volume from the total volume. The Thiessen polygon method is 
used to interpolate the rain gauge data to the spatial distributed grid data fitting the model 


123 grid resolution (30 arc-second) (Figure 1). We also obtained the daily discharge data at 


124 the basin outlet for the same time period. 

125 2.2.2 TMPA 3B42 Research Products 

126 TMPA 

’ (3B^ 


A precipitation products are available in two versions: near-real-time version 
127 (3B42RT) and post-real-time research version (3B42) adjusted by monthly rain gauge 


128 data. The 3B42 products have two successive versions: version 6 and the latest version 7 

129 (3B42V6 and 3B42V7). In this study, we evaluated and compared the two high-resolution 


130 (3 hours and 0.25°x0.25°) satellite precipitation products: 3B42V6 and 3B42V7. 
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131 The TMPA algorithm (Huffman et ah, 2007) calibrates and combines microwave (MW) 

132 precipitation estimates, and then creates the infrared precipitation (IR) estimates using the 

133 calibrated MW. After this, it combines the MW and IR estimates to create the TMPA 

134 precipitation estimates. MW data used in Version 6 are from the TRMM Microwave 

135 Imager (TMI), Special Sensor Microwave Imager (SSM/I) F13, F14 and F15 on Defense 

136 Meteorological Satellite Program (DMSP) satellites, and the Advanced Microwave 


137 Scanning Radiometer-Earth Observing System (AMSR-E) on Aqua, and the Advanced 

138 Microwave Sounding Unit-B (AMSU-B) N15, N16 and N17 on the NOAA satellite; IR 



139 data collected by geosynchronous earth orbit (GEO) satellites, GEO-IR. The 3B42V6 


140 also use other data sources: TRMM Combined Instrument (TCI) employed from TMI and 

141 PR, monthly rain gauge data from GPCP (l°xl°) and the Climate Assessment and 

142 Monitoring System (CAMS) 0.5° x 0.5 o developed by CPC. Based on the lessons learned 

143 in 3B42V6, 3B42V7 includes consistently reprocessed versions for the data sources used 

144 in 3B42V6 and introduces additional datasets, including the Special Sensor Microwave 

145 Imager/Sounder (SSMIS) FI 6- 17 and Microwave Humidity Sounder (MHS) (N18 and 

146 N19) and Meteorological Operational satellite programme (MetOp) and the 0.07° Grisat- 




147 B1 infrared data. All of these data can be freely downloaded from the website: 

X 

148 http://trmm.gsfc.nasa.gov/ and http ://mirador. gsfc .nasa. gov . 


149 2.2.3 Evapotranspiration 

150 The potential evapotranspiration (PET) data used in this study are from the global daily 


151 potential evapotranspiration database provided by the Famine Early Warning Systems 
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152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 


Network (hereafter referred as FEWSPET) global data portal (see 
http://carlvwarning.usgs.gov/fcws/global/wcb/readme.nhn7symboNDt) . FEWSPET is 
calculated from the climate parameter extracted from global data assimilation system 
(GDAS) analysis fields, has 1 -degree resolution, and covers the entire globe from 2001 to 
the present. 

2.3 CREST Model 

The Coupled Routing and Excess Storage (CREST) Model £ t ah, 201 la; Khan et 
ah, 2011b; Wang et ah, 2011) is a grid-based distributed hydrological model developed 
by the University of Oklahoma ( http://hydro.ou.edu) and NASA SERVIR Project Team 
( www.servir.net) . It computes the runoff generation components (e.g., surface runoff and 
infiltration) using the variable i n fi 1 trati ore ap acity curve (VIC), a concept originally 
contained in the Xinanjiang Model (Zhao, 1992; Zhao et ah, 1980) and later represented 
in the VIC Model (Liang et al., 1994; Liang et ah, 1996). Multi-linear reservoirs are used 



to simulate cell-to-cell routing of surface and subsurface runoff separately. The CREST 
model couples the runoff generation component and cell-to-cell routing scheme described 
above, to reproduce the interaction between surface and subsurface water flow processes. 
Besides the hydrologic and basic data (DEM, flow direction, flow accumulation, slope 
etc.), the CREST model employs gridded precipitation and potential evapotranspiration 
(PET) data as its forcing data. CREST Version 1.6 model has been applied at both global 
(Wu et al., 2012) and regional scales (Khan et al., 2011a; Khan et al., 2011b) (more 
applications can be found at website: http://eos.ou.edu and http://www.servir.net ). 
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173 The CREST model used in this study is the upgraded version CREST V2.0. The main 

174 features of the latest version are: 1) enhancement of the computation capability using 

175 parallel distribution techniques to make the model more efficient than the previous 

176 version (Wang et al., 2011); 2) model implementation with options of either spatially 

177 uniform, semi-distributed, or distributed parameter values; 3) automatic extraction of a- 

178 priori model parameter estimates from high-resolution land cover and soil texture data. 

179 The physically-based parameters, K sat and WM, can be derived from land cover types and 

180 soil texture data based on a look-up table (Chow et al., 1988); 4) a modular design 

181 framework to accommodate research, development and system enhancements (see Figure 

182 2 (a) ); and 5) inclusion of the optimization scheme SCE-UA (Duan et al., 1992; Duan et 


183 al., 1993) to enable automatic calibration of the CREST model parameters (see Figure 2 

184 (a) ). Table 1 shows 11 parameters and their descriptions, ranges and default values. 

185 Figure 2 (b) shows the vertical profile of hydrological processes in a grid cell. It shows 

186 the precipitation is intercepted by a canopy to generate throughfall, and then the 


6 


187 throughfall is separated into surface runoff and infiltration components by the variable 


188 infiltration curve. Finally, two linear reservoirs are employed to simulate sub-grid cell 

189 routing. 

190 


9 roul 

X 

o 


191 


Insert Figure 2 about Here 


192 


Insert Table 1 about Here 


193 
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194 2.4 Evaluation Statistics 

195 In order to quantitatively analyze the performance of 3B42V6 and 3B42V7 

196 precipitation products against rain gauge observations and the effect on 

197 simulation, three widely used validation statistical indices were selected in this 

198 relative Bias (%) was used to measure the agreement between the averaged value of 

199 simulated data (in this study, we call both TMPA products and simulated streamflow as 

200 “simulated data”, “SIM” was used in the formulae) and observed data (such as rain gauge 

201 and observed streamflow in this study, “OBS” was used in the formulae). The root mean 

202 square error (RMSE) was selected to evaluate the average error magnitude between 

203 simulated and observed data. We also use correlation coefficient (CC) to assess the 



204 agreement between simulated and observed data. 


205 


206 


r° 


o 




* 


Bias = 


'£sim,-'£obs, 

i=\ i = 1 

Z°BS, 


i = 1 


xlOO 


RMSE = 




i = 1 


n 

2 ( OBS . - OBS ) ( SIM, - SIM ) 
CC= i=l = 

^{oBs i -OBsf^ i [siM i -mi ) j 2 


( 1 ) 

( 2 ) 

( 3 ) 


V ;=1 i=l 

208 Where n is the total number of pairs of simulated and observed data; i is the rth values of 

209 the simulated and observed data; SIM and OBS are the mean values of simulated and 
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210 

211 

212 

213 

214 

215 

216 

217 

218 
219 


observed data respectively. Nash-Sutcliffe Coefficient of Efficiency (NSCE) is also used 
to assess the performance of model simulation and observation. 

V (OBS i -SIM.) 2 

NSCE = 1 - 

(« OBS-OBS ) 2 

3 Results and Discussion 

3.1 Comparison of Precipitation Inputs 

To better understand the impact of precipitation inputs on hydrologic models, the 
accuracy of the satellite precipitation against the in-situ rain gauge observations should be 
assessed first. This section compares the TMPA and gauge observations over the time 
span of January 1, 2001 to December 31, 2010 considering the basin-average 
precipitation and within a grid cell containing the dense rain gauge observations (Figure 





220 1). Figure 3 shows that both 3B42V6 and 3B42V7 systematically underestimate, though 


225 


226 


227 


B42V6 
'B42V7 

r 


221 at different levels, with biases of -41.15% and -8.38% and CCs of 0.36 and 0.40 at daily 

222 scale, respectively. Similar statistics are found at 0.25° grid-cell scale. Figure 4 indicates 

223 that 3B42V6 largely underestimates with a bias of -40.25% and low CC of 0.37, while 

224 3B42V7 has practically no bias (0.04%) and a relatively higher CC value 0.41. 


Insert Figure 3 about Here 


228 


Insert Figure 4 about Here 
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229 

230 Figures 3d and 4d present the inter-comparison of monthly precipitation estimates to 

231 gain further information about the precision and variations at longer time scales. The 

232 monthly data for both basin-based and grid cell-based analyses were accumulated from 

233 daily data over the same time span from January 2001 to December 2010. At monthly 


234 


235 


236 


237 


238 


239 


240 


time scale, both the basin-based and grid cell-based data show that 3B42V7 has better 

latest V7 algorithm significantly corrects the underestir bias in its predecessor 

version V6. 

Figure 5 (a) and (b) show the frequency distribution of daily precipitation for different 
precipitation intensities (PI) for the basin-averaged and the grid cell-based precipitation 
time series, respectively. Figure 5 (a) shows that for the basin- averaged data, both 



241 3B42V6 and 3B42V7 overestimate at the low PI range (less than 5 mm/day), but they 

242 underestimate at the medium and high PI ranges. However, 3B42V7 is in better 

243 agreement with the rain gauge observations than 3B42V6 for the basin-averaged 

244 comparison across all Pis. Similarly, better agreement has been found in Figure 5 (b) for 


245 the new Version-7 products at the grid cell scale, except for values greater than 30 

X 

246 mm/day where there is overestimation. 


247 


248 


Insert Figure 5 about Here 


249 
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250 3.2 Streamflow Simulation Scenarios 

251 Although different precipitation products vary in accuracy and spatiotemporal 

252 resolutions, they might have similar hydrological prediction (i.e., streamflow simulation) 

253 skill after re-calibrating the model using the respective precipitation products (Jiang et al., 

254 2012; Stisen and Sandholt, 2010). In the previous section, we compared the 3B42Y6 and 


255 


256 


257 


258 


259 


260 


261 


262 


264 


3B42V7 precipitation products against the rain gauge observations; the next step is to 


evaluate how these two TMPA products affect streamflow simulations. Their hydrological 
evaluation is performed under two scenarios: 

I. In-situ gauge benchmarking: Calibrate the CREST model with five years of rain 
gauge data (January 2001 through 



' simulation^ ' 
3REST n 


er 2005). Then, replace the rain gauge 
forcing with precipitation from 3B42V6 and 3B42V7 for an independent 
validation period from January 2006 through December 2010 using the rain 
gauge-calibrated model parameters. 


263 II. Product-specific calibration: Recalibrate the CREST model using the 3B42V6 and 


3B42Y7 precipitation data, respectively, over the same calibration period and then 


the product-specific parameter sets to simulate streamflow over the same 
validation period as Scenario I. 

267 Scenario I, gauge benchmarking, is widely used by the hydrological community 


' F Sc 


268 especially over gauged basins, while Scenario II is arguably deemed as an alternative for 

269 application to ungauged basins where only rainfall from remote-sensing platforms are 


270 available for use. 
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271 

272 3.2.1 Scenario I: CREST Benchmarked by In-situ Gauge Data 
1) Rain Gauge Calibration and Validation 

The CREST model parameters are calibrated using rain gauge inputs for the period 


273 


274 


A 


275 from January 2001 to December 2005 using the automatic calibration method (SCE-UA) 


276 by maximizing the NSCE value between the simulated and observed daily streamflow. 

277 The calibrated model is subsequently validated for the period from January 2006 to 

278 December 2010. Figure 6 compares the simulated streamflow forced by rain gauge data 


279 with the observed streamflow in terms of time series plots and exceedance probability 


; series 


280 plots at daily and monthly scales. Figure 6 (a) and (b) show that general agreement exists 

281 between the observed and simulated streamflow. However, the simulated streamflow 

282 consistently underestimates the peaks, especially in the validation period and in relatively 

283 low flow seasons as well. The exceedance probabilities in Figure 6 (c) and (d) also show 

284 underestimation at low and high streamflow observations, while the simulations match 

285 relatively well in the intermediate ranges. As summarized in Table 2 and Table 3, the 

286 statistical indices show that there is very good agreement between the simulated and 


287 observed hydrographs in the calibration period for both daily and monthly time scale, and 

r 

288 reasonable simulations occurred in the validation period as well. Based on the criteria of 


289 the statistical indices in Moriasi et al. (2007), the model calibration and validation results 

290 indicate that the CREST model is well benchmarked by the in-situ data at the daily and 

291 monthly time scale, so it can be used to evaluate the utility of the satellite precipitation 
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292 products for hydrological prediction (i.e., streamflow) in this basin. 


293 


294 


295 


296 


297 


298 


299 


Insert Figure 6 about Here 
Insert Table 2 about Here 
Insert Table 3 about Here 
2) Impacts of satellite precipitation forcing 

The gauge-benchmarked model is then forced by the TMPA 31 and 3B42V7 

A 

products from 2001 to 2010 using the model parameters by rain gauge data 

during the period from 2001 to 2005. Figure 7 and Fij iparc the daily and 



300 monthly time series of the simulated and observed hydrographs for both the calibration 

301 and validation periods. While 3B42V6 largely missed the high peak flows at both daily 

302 and monthly time series, 3B42V7 adequately captured a majority of the peak flows, 

303 especially at the smoothed monthly scale. The daily and monthly statistical comparisons 

304 in Table 2 and Table 3 show that the daily and monthly simulations forced by rain gauge 

305 data had better skill (NSCE=0.76/0.9 1 , BIAS=-9.73%/-9.75%, CC=0.89/0.96) than those 

306 based on 3B42V6 and 3B42V7 in the calibration period, which is expected. Interestingly, 


307 the 3B42V7-forced model simulations had very similar to and slightly better performance 

308 compared to the rain gauge-forced simulations in the validation period. A likely 

309 explanation is one of the rain gauge stations (i.e. the Dochula) had missing data from 


310 September 2006 to December 2010, which apparently degrades the hydrologic skill of 

311 this product. Overall, simulations forced by 3B42V7 were a significant improvement 


312 over 3B42V6. This clearly shows the improvements of the new version-7 algorithm upon 
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313 its predecessor V6 products both statistically and now hydrologically. 


314 


315 


316 


318 


Insert Figure 7 about Here 




* 


317 3.2.2 Scenario II: CREST calibrated by individual TMPA products 


X 


4 ? 

To further assess the effects of TMPA 3B42 (V6 and V7) products on streamflow, 


319 the CREST model is recalibrated and validated with 3B42V6 and 3B42V7 for the same 

320 period as Scenario I. This scenario is often used as an alternative strategy for remote 

321 sensing precipitation over ungauged basins. As shown in Figure 8, all simulations are 

322 significantly improved after the recalibration, and they capture most of the daily and 

323 monthly peak flows. Comparatively, the CREST model simulations based on 3B42V7 

324 inputs have better skill than those based on 3B42V6. As summarized in Table 2 and Table 


325 3, simulations have good statistical agreement with observed streamflow at daily and 

326 monthly scale. 

- 

328 visual comparison of the two modeling scenarios. Figure 9 indicates that the product- 

329 specific recalibration in Scenario II has obviously improved the NSCE and CC values 

330 and reduced the Bias values for both the calibration and validation periods. It is noted that 


331 the recalibration forcing with 3B42V7 in Scenario II has much higher NSCE and smaller 

332 Bias than 3B42V6, and very comparable CC values, all of which improved upon the rain 


333 gauge-benchmarked model. 
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334 


335 3.3 Discussion of parameter compensation effect from Scenario II 


336 


Table 4 shows the optimum parameter sets forced by 3B42V6 and 3B42V7, relative 


337 to the gauge forcing, for the calibration period from 2001 to 2005 using the SCE-UA 

338 algorithm. Note that the parameter values of Ksat and WM are spatially distributed but 

339 have been basin-averaged and summarized in Table 4. It shows that 3B42V7-calibrated 


340 parameters have less deviation from the gauge-calibrated para 


dues than 3B42V6. 


341 For example, RainFact is the adjustment factor of the 


ation either due to canopy 


342 


343 


344 


345 


interception or bias. Table 4 shows that 3B42V6 increases the RainFact value from 0.87 
to 1.34, to compensate its underestimation as shown in Figure 3 and Figure 4, while 
3B42V7’s estimated value (0.98) is closer to 1 and the Gauge value (0.87). Another 
example is KE, the ratio of potential evapotranspiration to the satellite PET data. Table 4 


346 reveals that 3B42V6 demands a reduced KE value from 0.10 to 0.05 in order to partition 

347 more precipitation into runoff while 3B42V7 only slightly increases it from 0.10 to 0.13, 

348 possibly to partially offset the above RainFact increase, amongst other parametric 


349 interactions. The third example is Ksat, the soil saturated hydraulic conductivity. Table 4 

350 shows that the Ksat of 3B42V6 reduced from 56.90 to 33.09 while V7 only changed 

351 slightly from 56.90 to 52.73. Regarding the mean water capacity, WM, 3B42V6 


352 decreased from 166.50 to 142.71 to hold less water in the soils while 3B42V7 did not 

353 change much from the gauge-calibrated value, which is presumably closer to the truth. It 

354 also shows the overland flow coefficient, coeM, the average channel flow speed, coeR, 
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355 the overland flow recession coefficient, KS, and the interflow recession coefficient, KI, 

356 all had reduced values to retain more water in the river basin after recalibrating the 

357 parameters to both of the satellite products. Not surprisingly, Table 4 also shows some 

358 opposite changes of values such as KE for 3B42V7 and coeS, the surface-interflow 


359 


360 


361 


362 


conversion factor, for both 3B42V6 and 3B42V7, resulting in a slight decrease in 

/V 

streamflow. 

In addition to the analysis of the parameters properties, water balance analysis is 
another important indicator for analyzing the effect of the parameter recalibration. Thus 


363 the difference of water balance components over 10-year (2001-2010) simulations is 

364 further examined using rain gauge and TMPA 3B42 rainfall, respectively. In CREST 

365 model, the water balance budgeting partitions the precipitation after canopy interception 

366 into actual evapotranspiration (ET), runoff depth (i.e. streamflow) and water storage 

367 change (AS), as shown in Figure 10. As expected, precipitation is the dominant runoff 

368 generation input so in Figure 3, all satellite rainfall forced simulations underestimated the 




369 streamflow compared to rain gauge results in scenario I. However, in scenario II, the 


370 model was recalibrated with respective satellite rainfall, the increased partition of the 
371, ^, s atellite driven streamflow simulations comes at the expense of a significant decrease of 
372 water storage due to the effect of the parameters value changes (shown in figure 10). In 


373 the gauge rainfall driven simulations, 27.90% of precipitation will be stored in this basin, 

374 however, 26.43% (27.18%) of precipitation is water storage in scenario I while 8.95% 


375 (16.09%) in scenario II for 3B42V6 (3B42V7). 
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From the above discussion, it is clear that the overall effect of the recalibrated 
parameter sets is to largely compensate for rainfall underestimation in 3B42V6 while less 
so for 3B42V7. The effect of arriving at a very similar simulation with different 


combinations of parameter settings has been called “Equifinality” of the hydrological 
model (Aronica et ah, 1998; Beven and Freer, 2001; Zak and Beven, 1999). This study 
clearly shows how different parameter settings can compensate for errors in the satellite 


rainfall forcing and can thus improve model predictions of : amflow. It is possible that 

the current model structural deficiency, i.e., not account! for snowmelt process, is 



compensated by the model re-calibration. However, this parameter compensation effect 
comes with the price of having a locally optimized model with parameter values 
unrepresentative of reality. This might limit the model’s predictive capability at internal 
sub-basins, or under different initial conditions. This is particularly concerning under 
scenarios involving climate change. In any case, the recalibration strategy could be 
especially problematic for 3B42V6 (Bitew and Gebremichael, 2011; Jiang et al., 2012), 
however the 3B42V7 product gives higher confidence for use in ungauged basins even 
without the need for recalibration. 


Insert Table 4 about Here 
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396 4 Summary and Conclusions 

397 Satellite precipitation products are very important for regional and global 

398 hydrological studies, particularly for remote regions and developing countries. This study 

399 first focuses on statistically assessing the accuracy of the TMPA 3B42V6 product vs. its 

400 latest successive version 3B42V7, and then hydrologically evaluates their streamflow 

401 prediction utility using the CREST distributed hydrologic model in the mountainous 


402 


403 


404 


405 


410 


411 


Wangchu Basin of Bhutan. 

The two versions of TMPA satellite products are statistically compared with a 
decade-long (2001-2010) rain gauge dataset at daily and monthly scales. In general, 
3B42V7 consistently improves upon 3B42V6’s underestimation both for the whole basin 



406 (bias improved from -41.15% to -8.38%) and for a 0.25° x 0.25° grid cell with high- 

407 density gauges (bias improved from -40.25% to 0.04%), though with modest 

408 enhancement of correlation coefficients (CC) (from 0.36 to 0.40 for entire basin and from 


409 0.37 to 0.41 fi 


or the 


lation ( 




grid cell). 3B42V7 also improves upon 3B42V6 in terms of 


occurrence trequency across the rain intensity spectrum. Apparently the results show that 
the new algorithm 3B42V7 has much improved accuracy upon 3B42V6, in concert with 


412 other studies in different areas (Chen at al. 2013ab and Kirstetter et al. 2013). The 

r 

413 improvement from V6 to V7 is mainly a combination of three factors: 1) the enhanced 


414 TMPA Level-2 retrieval algorithms (Chen et al. 2013a and Kirstetter et al. 2013), 2) 

415 incorporation of the global gauge network (i.e. GPCC) data with improved climatology 


416 and anomaly analysis (Huffman et al., 2011), and 3) additional satellite observations 
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417 incorporated (Huffman and Bolvin, 2012).. 


418 


421 


422 


423 


424 


425 


426 


427 


For the hydrological evaluation, two scenario-based calibration and validation 


419 experiments are conducted over the same 10-year time span. Scenario I, in-situ gauge 

420 benchmarking, is widely used by the hydrological community especially over gauged 


basins, while Scenario II, input-specific recalibration, is arguably deemed as an 
alternative for application to ungauged basins where only remote-sensing rainfall data 
may be available for use. In Scenario I, the 3B42V6-based si m ulation shows lower 
hydrologic prediction skill in terms of NSCE (0.23 at daily scale and 0.25 at monthly 
scale) while 3B42V7 performs fairly well (0.66 at daily scale and 0.77 at monthly scale), 
a comparable skill score with the simulations using the gauge benchmark. For the 
precipitation-specific calibration in Scenario II, significant improvements are observed 


428 for 3B42V6 across all statistics. These enhancements are not as obvious for the already- 

429 well-performing 3B42V7-calibrated model, except for some reduction in bias (from - 


430 26.98% to -4.81%). This behavior is consistent with previous studies (Bitew and 

431 Gebremichael, 2011; Bitew et al., 2011; Jiang et ah, 2012). This study offers unique 


432 insights into 3B42V6 and 3B42V7 products in a mountainous South Asian basin. 


43 W 

434 the l 


In concert with several other studies by Chen et al 2013a and Kirstetter et al 2013 in 
434 the US and Chen et al 2013b in the tropics, this study also reveals the latest 3B42V7 


435 algorithm has a noticeable improvements from 3B42V6 both in terms of accuracy (i.e., 

436 correcting the underestimation) and in its promising hydrological, even with or without 


437 recalibration of the hydrological model with respective rainfall inputs. The parameter 
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438 compensation effect is often recognized but still used by the hydrology community. This 

439 approach has been noted to be problematic due to unrealistic parameter settings which 


440 


441 


442 


446 


447 


448 


449 


450 


452 


453 


may ultimately limit the model’s predictive capability under conditions of climate change 
and differing initial conditions. 
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2 Figure 1 Wangchu Basin Map. (a) Location of Bhutan and the surrounding countries, (b) 

3 Location of Wangchu Basin in Bhutan and its elevation, (c) Map of Wangchu Basin, rain 

4 gauges, streamflow station, topography, Thiessen polygons applied to the rain gauge data 

5 and the 0.25°x0.25° grids of the satellite rainfall estimates. 
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13 Figure 3 Basin-averaged precipitation of a) Gauge, b) 3B42V6 and c) 3B42V7 for the 

14 period January 2001 -December 2010. Monthly data for both 3B42V6 and 3B42v7 are 

15 shown in d). 
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Figure 4 As in Figure 3, but for a single grid cell (Grid32). 
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Figure 5 Occurrence frequencies of rain gauge, 3B42V6 and 3B42V7 for a) basin- 

/rV 

averaged data and b) single grid cell (Grid32). 
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25 daily calibration (2001.1.1-2005.12.31) and validation (2006.1.1-2010.12.31), b) monthly 

26 data, c) exceedance probabilities using daily data from 2001.1.1 to 2010.12.31) and d) 

27 exceedance probabilities using monthly data. 
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Figur e 7 Comparison of CREST simulated streamflow from 3B42V6 (blue line) and 


B42V7 (red line) with gauge-calibrated parameters and observed stream flow in both 
calibration (2001.1.1-20.5.12.31) and validation (2006.1.1-2010.12.31) period, a) Daily 


32 data from 3B42V6; b) Daily data from 3B42V7; c) Monthly data from both 3B42V6 and 


33 


3B42V7. 
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Figure 8 As in Figure 9, but the parameters were recalibrated using 3B42V6 and 3B42V7, 
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respectively. 
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7 igure 9 Comparison of the streamflow performance statistics of the TMPA 3B42V6 and 
3B42V7 precipitation for the two simulation scenarios in both calibration period (a-d) 

and validation period (e-h). 
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Figure 10 Relative change of the water balance components using rain gauge and 


44 satellite rainfall based on ten-year annual averages (2001-2010) hydrological simulations 

cenarios I and II 
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1 Table 1 Monthly observed precipitation and runoff averaged from 2001 to 2010 


Rain Gauge (mm/month) 


Streamflow 

Station 


ivjL until 

Betikha Dochula 

Drukgyel Namjayling 

Dzong Haa 

DSC 

Paro 

Chhukha 
m 3 s' 1 

Jan 

14 

19 

0 

12 


8 

26 

Feb 

49 

18 

9 

24 


18 

fcb 

Mar 

Apr 

May 

176 

346 

368 

17 

53 

105 

26 

29 

60 

33 

54 

69 


15 

34 

57 

55 

Jun 

390 

279 

123 

124 


81 

111 

Jul 

546 

359 

185 

183 

i 

199 

222 

Aug 

383 

368 

191 

161 

103 

251 

Sep 

326 

217 

108 

12Qf\ 

77 

180 

Oct 

Nov 

Dec 

182 

10 

4 

116 

11 

9 

71 

3 

1 



63 

3 

1 

109 

52 

34 
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2 Table 2 Parameters to be calibrated in CREST V2.0, their description, ranges and default 


3 


values 


Parameter 

Ksat 

RainFact 

WM 

B 

IM 

KE 

coeM 

coeR 

coeS 

KS 

KI 


Description 


Numeric 

Range 

Default 

Value 

0-2827.2 500 

0. 5-1.2 1.0 

80-200 120 

0.05-1.5 T25 

0-0.2 

0.05 

0.1-1. 5 

1.0 

1-150 

90 

1-3 

2 

0.001-1 

0.3 

0-1 

0.6 

0-1 

0.25 


S 


soil saturated hydraulic conductivity (mm/d) 
multiplier on the precipitation field 
mean soil water capacity 
exponent of the variable infiltration curve 
impervious area ratio 

ratio of the PET to actual evapotranspiration 
overland runoff velocity coefficient 

multiplier used to convert overland flo 
speed to channel flow speed 

multiplier used to convert overland flow 
speed to interflow flow speed 

overland reservoir discharge parameter 
interflow reservoir discharge parameter 


o 
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5 Table 3 Comparison of daily observed and simulated streamflow under two calibration 

6 scenarios 


Precipitation Scenario I Scenario II 

Products NSCE Bias(%) CC RMSE NSCE Bias(%) CC RMSE 


Calibration Period 


Gauge 

0.76 

-9.73 

0.89 

45.38 

- 

- 


3B42V6 

0.23 

-52.94 

0.80 

81.99 

0.63 

-1.70 0.80 

56.55 

3B42V7 

0.66 

-26.98 

0.86 

54.65 

0.78 

-4.81 0.88 

43.94 




Validation Period 


o 


Gauge 

0.59 

-29.59 

0.83 

57.85 

- 

Ay 

- 

3B42V6 

0.17 

-57.78 

0.78 

82.98 

0.65 

-8.67 0.81 

54.00 

3B42V7 

0.63 

-25.15 

0.83 

55.26 

0.72 

-3.02 0.86 

47.72 





A 




r 



3 




ACCEPTED MANUSCRIPT 
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Table 4 As in Table 3, but for monthly data 


Precipitation Scenario I 


Scenario II 


Products NSCE Bias(%) CC RMSE NSCE Bias(%) CC RMSE 

Calibration Period 


Gauge 

0.91 

-9.75 

0.96 

25.18 

- 

- 

3B42V6 

0.25 

-53.01 

0.88 

72.08 

0.75 

-1.66 

3B42V7 

0.77 

-27.06 

0.94 39.76 

Validation Period 

0.91 

-4.83 

Gauge 

0.70 

-29.59 

0.88 

43.63 

- 


3B42V6 

0.19 

-57.81 

0.89 

71.35 

0.79 

-8.65 

3B42V7 

0.80 

-25.25 

0.94 

35.58 

0.89 

-3.08 


0.87 

0.95 


0.89 

0.95 




A 


36.29 

26.53 
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8 Table 5 CREST model parameter values calibrated with different precipitation inputs for 


9 


the calibration period of January 2001 -December 2005 


10 


Parameters Gauge 3B42V6 3B42V7 


RainFact 

0.87 

1.34 

0.98 

Ksat 

56.90 

33.09 

52.73 

WM 

166.50 

142.71 

166.52 

B 

1.48 

1.48 

1.48 

IM 

0.20 

0.20 

0.20 

KE 

0.10 

0.05 

0.13 

coeM 

88.05 

63.67 

67.95 

coeR 

coeS 

KS 

KI 

2.68 

0.43 

0.99 

0.20 

1.33 1.44 

0.47 0.67 

0.71 0.78 

0.13 0.14 


O 








5 




ACCEPTED MANUSCRIPT 


Highlights 


Comprehensively evaluated the latest satellite precipitation algorithm TMPA V6 and V7 


a 


A 


The new V7 provides more accurate spatiotemporal distribution than its predecessor V6 


V7 also provides improved hydrologic utility even without conventional recalibration 




ngaugeo I 

A 


& 


e$~ 


Promising potential application of the latest V7 in ungauged basins across the world 





